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ABSTRACT 



The upcoming generation of SZE surveys will shed fresh light onto the 
study of clusters. What will this new observational window reveal about cluster 
properties? What can we learn from combining X-ray, SZE, and optical observa- 
tions? How do variations in the gas entropy profile, dark matter concentration, 
accretion pressure, and intracluster medium (ICM) mass fraction affect SZE 
observables? We investigate the signature of these important cluster parameters 
with an analytic model of the ICM. Given the current uncertainties in ICM 
physics, our approach is to span the range of plausible models motivated by 
observations and a small set of assumptions. We find a tight relation between 
the central Compton parameter and the X-ray luminosity outside the cluster 
core, suggesting that these observables carry the same information about the 
ICM. The total SZE luminosity is proportional to the thermal energy of the 
gas, and is a surprisingly robust indicator of cluster mass: L$z oc /icm^ 5//3 - We 
show that a combination of L$z and the half-luminosity radius rsz provides a 
measure of the potential energy of the cluster gas, and thus we can deduce the 
total energy content of the ICM. We caution that any systematic variation of 
the ICM mass fraction will distort the expected L sz — M calibration to be used 
to study the evolution of cluster number density, and propose a technique using 
kSZ to constrain ficM^M, z). 
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Subject headings: cosmic microwave background — cosmological parameters - 
galaxies: clusters: general — X-rays: galaxies: clusters 

1. Introduction 

As the largest and most recently formed relaxed objects in the universe, clusters pro- 
vide cosmological constraints complementary to other observations. Press & Schechter (1974) 
showed that their co-moving number density is exponentially sensitive to both cluster mass 
and the amplitude of the linear power spectrum of density fluctuations. The number density 
of massive galaxy clusters in the local universe constrains a model-dependent combination 
of the amplitude of fluctuations at 8 h" 1 Mpc (<j 8 ) and the matter density (f2 m ). The red- 
shift evolution of cluster number counts depends on the linear growth factor and and the 
co-moving volume element, and can thereby potentially break the degeneracy in the family 
of A and quintessence models allowed by primary CMB anisotropy (Wang & Steinhardt 
1998). X-ray and optical cluster counts at low and intermediate redshifts have been used 
to constrain cosmological parameters (for recent results, see Allen et al. 2003, 2004; Bahcall 
et al. 2003; Henry 2004; Ikebe et al. 2002; Pierpaoli et al. 2003; Rapetti et al. 2005; Reiprich 
& Bohringer 2002; Rosati et al. 2002; Shimizu et al. 2003; Viana et al. 2002; Vikhlinin et al. 
2003). 

The Sunyaev-Zel'dovich effect (SZE) flux is an excellent probe of high redshift clusters, and 
many SZE surveys are in development, such as ACT 1 , APEX 2 , Planck 3 , and SPT 4 . While 
the information from high redshift available in SZE surveys can potentially provide tight 
cosmological constraints (Haiman et al. 2001), recent work suggests that uncertainties in the 
mass-observable relation and its scatter can severely degrade their sensitivity. Fortunately, 
internal calibration and some follow-up observations can recover much of a survey's sensitiv- 
ity if the mass-observable relation and its scatter can be described accurately in the relevant 
redshift range by a reasonably small set of parameters (Hu 2003; Majumdar & Mohr 2004; 
Lima & Hu 2004, 2005; but see also Francis et al. 2005). The results presented in this paper 
support their assumption: we find a tight mass-observable SZE relation that is robust to a 
wide range of model variations. 

The SZE is a distinct spectral signature in the CMB due to Thomson scattering of CMB 

1 http: / /www. hep. upenn.edu/act / 
2 http://bolo. berkeley.edu/apexsz/index. html 

3 http: / /www.rssd.esa.int /index. php?projcct==PLANCK&page=indcx 
4 http://spt. uchicago.edu/ 
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photons and hot electrons with a magnitude proportional to the integrated electron pressure 
along the line of sight (for discussions of SZE physics, see Birkinshaw 1999; Carlstrom et al. 
2002). A unitless measure of the effect is the Compton parameter y (Sunyaev & Zeldovich 
1972): 

^Thomson / n i~ \ 

y = — / Pe dl. 1 

m e c 2 J 

We define the SZE luminosity Lsz as the integration of y over the projected cluster area, so 
that 

/sz4 = L SZ = [ydA= f P e dV = ^^pN e (T e ). (2) 

J m e c 2 J m e c 2 

Assuming local thermal equilibrium, L$z directly measures the total thermal energy content 
of the intracluster medium (ICM), and y(8) informs us on how the gas is distributed in the 
cluster. Because cLa flattens at intermediate redshift, the SZE flux / sz becomes approxi- 
mately redshift independent. 

The self-similar collapse model for the ICM considers only gravitational physics and predicts 
scaling relations between cluster properties. X-ray observations indicate significant devi- 
ations from these predictions, most notably in the X-ray luminosity-temperature relation 
(Markevitch 1998; Arnaud & Evrard 1999). Many complex non-gravitational processes de- 
termine the final thermal state of the ICM. We discuss many of them in Sec. 2 along with 
related X-ray observations. In Sec. 3 we introduce an analytic model of smooth accretion 
(Voit et al. 2003) that provides a conceptual foundation for determining the gas properties. 
We then relax some of the assumptions of Voit et al. (2003) and parametrize our uncertain- 
ties in ICM physics: the poorly understood effects of heating, cooling, and conduction is 
encoded in the gas entropy profile and the ICM mass fraction, variation in the dark matter 
potential is characterized by the concentration parameter, and the boundary condition is 
set by an accretion pressure. Motivated by the physics and observations discussed in Sec. 2, 
our main goal is to explore SZE and X-ray observables over a plausible range of parameter 
values using a set of phenomenological models under the assumptions of spherical symmetry 
and hydrostatic equilibrium. We compare the resulting profiles to recent observations of 
nearby, relaxed clusters reported in Vikhlinin et al. (2005). We quantify deviations from the 
expected L sz — M relationship and indicate the main sources of scatter accessible to our 
models. In Sec. 4 we describe relations between cluster properties and observables that hold 
across the range of models we consider. In particular, we show that SZE observables alone 
provide a robust measure of the total energy content of the ICM. We discuss the implications 
and limitations of our models in Sec. 5, and focus on the important point of constraining the 
ICM cluster mass fraction /icm using the kinetic Sunyaev-Zel'dovich effect (kSZ) in Sec. 6. 
Sec. 7 states our conclusions and the implications of this work for future SZE surveys. 
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2. ICM Physics 

In the self-similar spherical collapse model (Peebles 1993) and in agreement with N-body 
simulations, all virialized clusters at fixed redshift have the same overdensity relative to the 
background. In the absence of non-gravitational processes, we expect the gas temperature to 
be set by the dark matter virial temperature: T vir = GM ^ n f qq M 2//3 where fj,m p is the mean 
gas particle mass. These assumptions also predict Lx oc T 2 when thermal bremsstrahlung 
dominates X-ray emission (T > 2 keV) and L x oc T at lower temperatures when line emis- 
sion dominates. However, observations indicate Lx oc T 2 - 6-2 - 9 (Markevitch 1998; Arnaud 
& Evrard 1999). This deviation from the expected scaling has instigated extensive theo- 
retical investigation of the effects of various non-gravitational processes. Preheating of the 
intergalactic medium at an early time (Kaiser 1991), radiative cooling (Voit & Bryan 2001), 
cooling with energy injection due to feedback (Ostriker et al. 2005), and quasar activity 
(Lapi et al. 2005) can all reproduce the observed L x — T relation. Fortunately, the latest 
high resolution X-ray observations provide new constraints for ICM models out to nearly 
half the dark matter virial radius (Vikhlinin et al. 2005). 

Under the assumptions of spherical symmetry and hydrostatic equilibrium, X-ray surface 
brightness and temperature profiles allow extraction of a cluster's total mass profile. Recent 
data from nearby, relaxed clusters (Vikhlinin et al. 2005; Pratt & Arnaud 2005) are consis- 
tent with the universal NFW profile (Navarro et al. 1997) with ACDM best fit concentration 
parameters and scatter (Dolag et al. 2004). In Vikhlinin et al.'s sample of 13 clusters with 
temperature range 0.7-9 keV, temperature profiles are self-similar only at r > 0.15r v i r . 
However, using the gas mass weighted rather than X-ray emission weighted temperature, 
they still find agreement with self-similar predictions for the M — T scaling. 
Since cluster gas evolves nearly adiabatically, the entropy profile is a useful description of 
the gas: once the entropy distribution is fixed, convective stability, a confining pressure at 
the cluster gas boundary, and the assumption of hydrostatic equilibrium determine the gas 
profile in a dark matter potential. The customary definition of 'entropy' in this field is (Tozzi 
& Norman 2001; Voit et al. 2003) 

K = = T oc e 2s/3 (3) 
p 5 / 3 firripp 2 / 3 

and is related to the thermodynamic entropy per particle s as given above. Typically, one 
considers the function K(M gas ), where M gas is the mass of gas with entropy is < K(M gas ). 
A useful quantity is the characteristic entropy of a region of overdensity A as defined by 



(4) 
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with Ta = — 2 ^ mp the characteristic temperature of the dark matter potential at overdensity 
A, and /b the universal baryon fraction. The self-similar collapse model predicts a universal 
entropy profile K(r) oc r 11 fixed by the physics of gravitational collapse (Tozzi & Norman 
2001; Borgani et al. 2001), with the normalization scaling as K\ oc Ta for clusters at fixed 
redshift. Observed entropy profiles (Pratt & Arnaud 2005) are approximately self-similar 
down to 2 keV in the radial range 0.05r 2 oo < r < 0.5r 2 oo with K(r) oc r°- 94±014 (Pratt & 
Arnaud 2005), and they generally show a flatter core inside 0.1r 2 oo (Pratt & Arnaud 2003; 
Ponman et al. 2003). The entropy profile normalization clearly deviates from the self-similar 
prediction and scales with cluster temperature as i^(0.1r 20 o) oc T a65±a05 (Ponman et al. 

2003) . 

As suggested by Neumann & Arnaud (2001), similarity breaking in both the Lx — T and 
K — T relations are consistent with an ICM mass fraction temperature dependence, so that 
/icm ~ T 1 / 2 . X-ray observations find evidence for lower gas densities in cooler clusters 
out to ~ 0.35r 20 o (Neumann & Arnaud 2001; Pratt & Arnaud 2003; Vikhlinin et al. 2005; 
but see also Mathews et al. 2005). Current X-ray observations cannot measure gas profiles 
much past 0.5r vir , so the 'missing' gas may be bound in the outer regions of the cluster, 
escaped from the cluster potential altogether, or condensed into intracluster stars (Gonzalez 
et al. 2005; Lin & Mohr 2004). SZE signatures in WMAP also find that /icm increases with 
temperature (Afshordi et al. 2005). Because L sz oc /icm, any trend in the total ICM mass 
fraction will affect the Lsz — M relation. Such trends are quite possible, given trends with 
mass in observations of the cluster mass/near- infrared K-band luminosity relation (Lin et al. 

2004) , estimates of intracluster light (Lin & Mohr 2004), and in hydrodynamic simulations 
finding a mass-dependent ICM mass fraction (Kravtsov et al. 2005). 

Voit & Bryan (2001) have suggested that radiative cooling sets the entropy scale responsible 
for similarity breaking: gas below the cooling threshold is either condensed or injected with 
energy through feedback until it exceeds the cooling threshold. Consistent with observed 
scaling relations, their simple model predicts for the central entropy K oc T 2//3 for T > 2 keV 
and K independent of T for T < 2 keV. Dave et al. (2002) demonstrate a similar K — T 
scaling in numerical simulations of groups that include radiative cooling and star formation. 
In the analytic model of McCarthy et al. (2004) cooling approximately maintains the initial 
entropy power law, in agreement with the observed entropy gradients in clusters. 
Cooling and subsequent feedback play an essential but enigmatic role in the thermodynamics 
of the ICM. As many as 70 — 90% of nearby clusters have cold cores of gas with t coo \ much 
smaller than the age of the cluster (Peres et al. 1998), and similar results have been found in 
the redshift range z xs 0.15 — 0.4 (Bauer et al. 2005). However, in several well-studied cold 
core clusters, Tamura et al. (2001), Peterson et al. (2001), Sakelliou et al. (2002), Peterson 
et al. (2003), and Kaastra et al. (2004) find very little cooled gas below a quarter of the hot 
gas temperature. Measurements of O VI emission in cold cores imply gas condensation rates 
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well below the expectations from a simple cooling flow model (Bregman et al. 2006), and 
altogether observations suggest subsonic energy injection spatially distributed in regions of 
~ 100 kpc and continuous on timescales of 10 8 years. Thermal conduction and AGN activity 
may be important to balancing radiative cooling in cold cores (see Peterson & Fabian 2005, 
for a comprehensive review of the cooling flow problem). 

A distinct source of feedback is needed to balance the severe overcooling present in cooling- 
only analytic calculations and cosmological hydrodynamic simulations (Oh & Benson 2003; 
Balogh et al. 2001). The stellar baryon fraction in clusters has been estimated at 5 — 20% 
(Roussel et al. 2000; Balogh et al. 2001; Lin et al. 2003), though inclusion of intracluster 
stars may increase this fraction by a factor of two (Lin & Mohr 2004; Gonzalez et al. 2005; 
Zibetti et al. 2005). Voit & Bryan (2005) conclude that heating is needed at high redshift to 
solve this 'cooling crisis,' and its resolution is most likely related to the high entropy levels 
observed in less massive clusters and groups. 

A simple preheating of the intergalactic medium to a 'universal entropy floor' by supernovae 
and AGN before the epoch of cluster formation (Kaiser 1991) introduces deviations from the 
expected scaling relations and mitigates the cooling crisis. Heating occurring before accre- 
tion onto the cluster most efficiently raises the gas entropy level. Lapi et al. (2005) estimate 
1/4 keV/particle is available from supernovae, and is not enough to match the observed 
Lx — T and K — T relations. Lapi et al. (2005) also consider an additional 1/2 keV/particle 
available for preheating from smooth, long-lived AGN outputs; this energy is sufficient to 
provide marginal agreement with the observed Lx~T and K — T scaling relations. However, 
external preheating models with smooth accretion are generically disfavored observationally 
because they predict large isentropic cores in low mass systems (Tozzi et al. 2000; Pratt & 
Arnaud 2003; Ponman et al. 2003) and require a large entropy floor (> 3x 10 33 erg cm 2 g~ 5 / 3 ) 
to match the L x — T relation (McCarthy et al. 2004). We do consider this model in our 
computations for completeness. 

Simulations show that the accretion process is far from spherical, often taking place along 
filaments. Thus many researchers (Voit et al. 2003; Ponman et al. 2003; Voit & Bryan 2005) 
have suggested a more favorable scenario of preheating where the effect of preheating is pri- 
marily to reduce the density of pre-shock gas accreting onto lower mass systems. This raises 
the post-shock entropy normalization while producing the radial entropy profile expected 
from gravitational shock heating. Borgani et al. (2005) have explored this proposal in nu- 
merical simulations, and find that feedback must be better distributed than their strongest 
galactic winds in order to produce the desired levels of entropy amplification and radial gradi- 
ent. Lapi et al. (2005) circumvent isentropic cores and introduce similarity-breaking through 
quasar blast waves within the cluster, which leave a steep final entropy profile K(r) oc r 1 ' 3 . 
In Sec. 3 we incorporate features from all the models discussed above into our set of phe- 
nomenological models. 
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3. SZE Signals from Analytic Models 

In this section we first describe our analytic ICM model taken from Voit et al. (2003) 
with slight modifications. Under the assumption of smooth accretion and gravitational shock 
heating, the mass accretion history determines K(M g ), which can be modified to account 
for cooling and a uniform preheating. Our phenomenological approach generalizes this ICM 
model to explore parameter ranges suggested or unconstrained by observations. We no longer 
assume smooth accretion or attempt to compute K(M g ) or /icm from first principles, and 
we also allow the dark matter potential and boundary pressure to vary. Because researchers 
(e.g., Ostriker et al. 2005) sometimes use a polytropic model for the ICM, we investigate its 
SZE observables and use it to consider a model of lumpy accretion. Throughout these models 
we translate the effects of non-gravitational physics and other uncertainties in ICM physics 
into their effects on the gas entropy profile, cluster potential, and boundary conditions, from 
which we can easily compute the cluster observables. Following Voit et al. (2003), we assume 
for cosmological parameters t Q = 13.4 Gyr, fi mj0 = 0.3, VL\ yQ = 0.7, h = 0.71, a 8 = 0.9, and 
the universal baryon fraction /b = 0.02/i~ 2 fi~ 1 o . We do not expect these choices to limit the 
validity of our results. In particular, the universal baryon fraction only sets the normaliza- 
tion of the gas density, since we neglect its self-gravitational effects. We also assume a fully 
ionized ICM with primordial abundances, so p/nm p = ji = 0.59 and n e = 0.52n. For cooling 
computations we assume a metallicity Z = 0.3Z Q . 



3.1. Analytic Model of Preheating and Cooling 

We model the dark matter with an NFW potential (Navarro et al. 1997) of overdensity 
A(z), mass M A , and virial radius r A : 

GM A log(l + c A r/r A ) r A ,_ N 

<PNFw{r) = : -— 7 7-— r • (5) 

r A log(l + c A ) - c A /(l + c A ) r 

We fix A(z) using the approximation in Bryan & Norman (1998) to the spherical top hat 
collapse model in a AC DM universe, so A ^ 100 at z — 0. As in Voit et al. (2003), we 
parametrize the concentration parameter's weak dependence on cluster mass and redshift 
(Navarro et al. 1997) as 

CA (M A ) = 8.5 (Ma/IO 15 /*-^)- - 086 (1 + z)~ - 65 , (6) 
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and we use the mass accretion histories fit in Voit et al. (2003) for clusters observed at z = 
to a merger tree algorithm of Lacey & Cole (1993): 

log(M/M D ) = ai log(f/fo) + a 2 (log(t/t )) 2 . (7) 

These functions already introduce some similarity breaking in our baseline model because a±, 
a,2, and ca all depend on cluster mass. We consider 1O 15 /i _1 M (ai = 1.94 and a 2 = —0.55), 
1O 14 /i" 1 M ( ai = 1.10 and a 2 = -0.88), and 1O 13 /i _1 M (a x = 0.64 and a 2 = -0.96) clusters 
at z = (Voit et al. 2003). 

In the smooth accretion approximation, infalling gas of uniform density is shocked at a 
well-defined radius, inside of which we assume hydrostatic equilibrium. The parameter £ 
describes the position of the gas accretion shock radius r ac 

e = i-^. (8) 

7"ta 

We assume r ta = 2r A (Voit et al. 2003) for the position of the turn-around radius of the 
gas currently accreting. We follow Voit et al. (2003) for a self-consistent determination of 
£(t) 5 . As in Tozzi & Norman (2001) and Voit et al. (2003), £ is nearly constant in time, 
except for a sharp jump as accretion goes from adiabatic to shock-dominated in preheating 
models. The free-fall velocity of the gas relative to the cluster at the accretion shock is fixed 
by energy conservation: 

<4 - « (9) 

' ac 



Note that we neglect gas energy changes associated with adiabatic compression of the in- 
falling gas and a time-varying amount of enclosed dark matter (see Tozzi & Norman (2001) 
for a discussion of these effects) . We adopt the approximation that dark matter and baryonic 
mass densities trace one another before accretion: 

f h M'(t) = 47rr a 2 c pi«i. (10) 

Here /b is the universal baryon mass fraction, r ac = 2rA(l — £) is the gas shock radius, p\ 
is the average pre-shock gas density, and u\ is the velocity of the gas relative to the moving 
shock radius. As our self-consistent determination of £ suggests, we further assume that the 
ratio of the shock radius to the dark matter virial radius changes slowly, so that U\ is just 
given by 

Ul= Vac + 2(1-0^(0- (11) 

Pi, Ui, and the external gas entropy i^ prc heat fix the Mach number of the shock at r ac . We 
apply the Rankine-Hugonoit jump conditions (Landau & Lifshitz 2000) to determine the 



5 We allow their parameter gi to vary and also approximate r' A (t) = 0.2v ac (t). 
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post-shock density and pressure. Thus M(t) and £(t) fix the post-shock entropy profile 
K(M gas ). Note that we do not assume that the post-shock velocity is zero, since the gas 
is continually compressing to maintain an approximately fixed value of £. The post-shock 
velocity relative to the cluster is < 25% of the sound speed throughout the accretion histories, 
so that hydrostatic equilibrium remains a good approximation. Our model ignores any 
subsequent accretion shocks bringing the gas to rest in the cluster; these would occur at 
higher densities and involve relatively small changes in energy. 

Once we have the entropy distribution K(M gas ) that includes preheating and accretion shock 
heating, we modify K(M gas ) according to the cooling approximation in Voit et al. (2003). 
Finally we solve the equations of hydrostatic equilibrium and mass conservation to determine 
the gas property profiles 



ldP d<t> 



NFW 



p dr dr 

= Airr 2 p 12 

dr 

P = K(M gas )p 5 / 3 . 



We vary £(t ) until M gas (r ac ) = ficMM(t ), where the mass accretion rate, the shock Mach 
number, and K(M(t )) fix the boundary conditions. To avoid the singularity associated 
with K = when the gas cools, we set a minimum entropy K min = O.Ol^oo so that we do 
not artificially introduce similarity breaking. SZE cluster properties are unaffected by this 
choice. 

To relate our gas profiles to an X-ray observable, we choose to compute Lx, C ut, the X-ray 
luminosity outside a projected radius 0.05rA, to avoid uncertainties about the small entropy 
values in the model cluster cores. We also report the total X-ray luminosity, L x , and the 
radius projected onto the sky containing half the total X-ray luminosity, rx (including the 
region inside 0.05t-a). The X-ray luminosity of a small region of gas dV is given by 

dL x = n 2 (x)A(T)dV. (13) 

A(T) is the X-ray cooling function, which we approximate with the cooling function that 
includes all wavelengths given in Sutherland & Dopita (1993). We compute SZE profiles y(9) 
as well as the SZE luminosity Lsz, the central Compton parameter y (y along the line of 
sight passing through the cluster center), and the projected radius containing half the SZE 
luminosity, r sz . 
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3.2. Preheating and Cooling Model Results 

We first compute the predictions of the model described in Sec. 3.1 when both preheat- 
ing and cooling are neglected, and then include cooling and allow i^ pre h e at to vary in the 
range 10 30 — 10 34 erg cm 2 g~ 5 / 3 . In Appendix A we report X-ray and SZE observable prop- 
erties of our model clusters. Our baseline model neglects cooling and sets i^ pre heat = 0, and 
produces an X-ray core surface brightness too large compared to observed cold core regions 
(Markevitch 1998). The systematic decrease in ca and increase in accretion pressure with 
increasing cluster mass already introduce some similarity breaking: low mass clusters are 
hotter in the cluster interior and have larger accretion radii. Otherwise gas profiles remain 
approximately self-similar. Our baseline model reproduces the self-similar Lsz scaling result, 
and we find L sz = (1.6 x 10" 4 Mpc 2 )(M/lO 15 /i^ 1 M ) 5 / 3 . 

With cooling but negligible pre-shock gas entropy, the cooled fraction of baryons was 0.05, 
0.15, and 0.42 for the highest to lowest mass clusters. The cooling approximation of Voit 
et al. (2003) depends on the assumption of T = T A throughout the cluster at all times, 
and appears to underestimate the cooled fraction compared to more precise analytic cooling 
models by around a factor of 2 (see Oh & Benson 2003, Fig. 1), and more so in cooling-only 
hydrodynamic simulations (da Silva et al. 2004). Cooling steepens the inner entropy profile 
sufficiently to produce X-ray core surface brightness values in the range measured in the 
cold core regions of Markevitch (1998). Once scaled according to the remaining gas fraction 
/icm, density and entropy profiles remained approximately self-similar. The Lsz — M scaling 
relation slightly steepens to L sz = (1.6 x 10~ 4 Mpc 2 )(M/lO 15 /i~ 1 M ) L75 . If we scale L sz 
by the uncooled gas fraction /icm as well, we recover the self-similar result found in the no 
cooling model above: L sz /(f 1CM / f h ) = (1.6 x 10~ 4 Mpc 2 )(M/lO 15 /i- 1 M ) 5 / 3 . Thus, in our 
model the mass-weighted temperature stays fixed when we allow cooling. 
In addition to cooling we now consider -Kp re heat in the range 10 30 — 10 34 erg cm 2 g~ 5 / 3 . Voit 
et al. (2003) adopt 3 x 10 33 erg cm 2 g -5 / 3 to reproduce the Lx — T relation, and obser- 
vations suggest 1.4 x 10 33 erg cm 2 g~ 5 / 3 (Lloyd-Davies et al. 2000; Ponman et al. 1999). 
Large Lf pre heat values introduce isentropic, hot cores, particularly in the least massive clus- 
ter. Figure 1 compares the SZE surface brightness profiles produced by the cooling-only and 
-^preheat = 3 x 10 33 erg cm 2 g~ 5 / 3 models. Preheating dramatically broadens SZE profiles 
in the low mass systems due to the fractionally greater energy input, while the Lsz scaling 
only slightly steepens for this large K prc hcat value: L sz = 1.6 x 10~ 4 Mpc 2 (M/lO 15 /i- 1 M ) L7 °. 
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Fig. 1. — SZE surface brightness profiles for our cooling only and cooling+preheating 

(-Kprcheat = 3 X 10 erg CHI g 

-5/3) mo( jel s scaled by M/lO 15 /i~ 1 M ( see Sec. 3.2 for more 
discussion). In a completely self-similar model the profiles of all cluster masses would coin- 
cide. Our cooling model (thin curves) breaks similarity by inducing a mass dependent cooled 
gas fraction, while our cooling+preheating model (thick curves) produces more diffuse cores 
in lower mass clusters. The solid curves are the 1O 15 /i _1 M cluster models. The primary 
result of preheating for this cluster is to shut off cooling, so the preheating model has larger 
gas fraction and thus SZE profile. Short-dashed curves are the 1O 14 /i _1 M cluster models, 
where preheating has significantly reduced the central gas density and therefore y a . The 
long-dashed curves are the 10 13 /i _1 M Q cluster models. The large isentropic core severely 
extends the gas profile in the preheating model. 
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3.3. Phenomenological Models and Comparison with Observations 

In this section we adopt a phenomenological approach in order to include information 
from the latest X-ray observations and to explore a wider range of cluster parameters. While 
we still assume spherical symmetry, an NFW potential (Eqn. 5), hydrostatic equilibrium, 
and a well-defined accretion shock radius bounding the ICM, we do not attempt to compute 
K(M gas ) or /icm from first principles. Furthermore, we allow the accretion pressure and 
dark matter concentration to vary. Motivated by both observations and models discussed 
in Sec. 2, we parametrize the gas entropy by a normalization K max and profile f(x = r/r&) 
with core radius x c and two exponents Si and s 2 ; 



A power-law parametrization is consistent with both the predictions from gravitational heat- 
ing and quasar blasts and with X-ray observations of entropy profiles (see Sec. 2). A flatter 
entropy profile (s 1 < s 2 ) in the core is consistent with observations (Pratt & Arnaud 2003; 
Ponman et al. 2003) and preheating models. 

Solving the hydrostatic equilibrium equation for the normalized pressure profile p(x) = 
P(x)/P(r ac ) with x = r/r& we find 



with T max = /j,m p P(r ac ) 2 / 5 K?neL, the gas temperature at r ac . We fix P(r ac ) by the gas ac- 
cretion rate ficMM'(t), incoming velocity U\ (Eqn. 11), and assuming a strong shock at t Q 
(that is, the external thermal pressure can be neglected at t Q ). We vary r ac until the en- 
closed gas mass is M gas = /icm^a- We do not assume smooth accretion, so if max varies 
independently of P(r ac ). Note that for fixed T max , the pressure, density, and SZE profiles 
are fixed, and their normalizations simply scale with /icm- This scaling will be broken 
only if cluster physics condenses or ejects a different fraction of gas compared with the 
local infalling average. A fit to the accretion rates at t from Voit et al. (2003) yields 
M'(t ) = (1.94 x lO^h^MQt-^dMiQ/lO^h^MQf^w^, where we allow a "fudge fac- 
tor" w accr in the accretion rate. 

In our "self-similar" model we set K max = K 100 , K min = 0.01K WQ , si = s 2 = 1.1, x c = 0, ca 
according to Eqn. 6, and /i CM = 0.13. These are the expected values in the gravitational- 
only heating scenario described in Sec. 2 (recall A(z = 0) = 100 in AC DM), with further 
cooling and condensation of baryons in clusters neglected. Parameters in other models are 
fixed to these values unless specified otherwise. We list the computed observables from our 
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models in Appendix A. Figure 2 and Figure 3 show the density and temperature profiles for 
phenomenological models of the most massive cluster; the others are similar. 
According to AC DM numerical simulations (Dolag et al. 2004), dark matter concentration 
parameters show intrinsic scatter in a log-normal distribution with ln(c/c avg ) ~ 0.22, and a 
range of mean values have been reported (Navarro et al. 1997; Komatsu & Seljak 2001; Dolag 
et al. 2004). We allow ca to range by 8 to encompass these uncertainties, and find a ~ 20% 
decrease in L sz for the lowest values of ca and « 10% increase for the highest for all cluster 
masses considered. Taking the mean ca and its scatter reported by Dolag et al. (2004) and 
our self-similar gas model, we find that the intrinsic scatter in ca induces variations in Lsz 
of < 8%. The effect is smaller at lower concentrations (higher mass clusters) and shallower 
entropy profiles (lower si). Figure 2 and Figure 3 show that increasing ca both concentrates 
more gas in the inner, hotter cluster region and significantly increases the temperature in 
that region. 

Mass accretion is stochastic, aspherical, and has not been directly observed at the cluster 
boundary. We vary the accretion pressure by a factor u> aC cr = 3.5 and w accr = 1/3.5. The 
additional accretion pressure boosts L S z by 44% in our largest mass cluster and pushes the 
gas temperature well above T A at r ac (see Figure 3). All other cases saw much milder changes 
in L S z- For the lQ lA h~ l M & and 1O 13 /i _1 M clusters we also list the results that would be 
obtained in the case of a self-similar mass accretion rate (M'(t ) oc M(t )). 
We allow for a core in the entropy profile, as suggested by preheating models and observed 
profiles, but assume a power-law entropy profile outside the core. The entropy profile is 
observationally unconstrained near the cluster virial radius, but gravitational heating prob- 
ably dominates in that region. We allow the power law index to range from 0.7 to 1.5, 
incorporating the observed K(r) oc r a94±a14 (Pratt & Arnaud 2005), K(r) oc r 11 ex- 
pected from gravitational heating (Tozzi & Norman 2001; Borgani et al. 2001), and the 
steep profile K(r) oc r 13 expected from quasar blasts (Lapi et al. 2005). We vary K max 
over more than a factor of 3, which should bound the expected range for clusters be- 
tween 1O 13 /i _1 M and lO 15 /i _1 M . Variations on S2, inclusion of an entropy core, fixing 
-f^min = 3.5 x 10 33 erg cm 2 g _5,/3 , and varying K max between i^oo and 3i^ioo all caused < 5% 
variation of Lsz for our M = 1O 15 /i _1 M cluster. Even though increasing K mSiX raises the 
temperature profiles, the gas also becomes more extended and occupies the cooler regions of 
the cluster. Our 10 14 /i _1 M o and 1O 13 /i _1 M clusters were slightly more sensitive to changes 
in the entropy normalization, but Lsz still remained within « 10% of the self-similar model. 
We consider two scenarios consistent with the observationally suggested K ma , x oc T 2 / 3 (Pon- 
man et al. 2003). If the primary effect of preheating is to reduce the pre-shock gas density 
by a factor dependent on cluster potential depth (Voit et al. 2003; Ponman et al. 2003; Voit 
& Bryan 2005) or if the fractional feedback entropy injection into the bound gas depends on 
cluster mass, then we expect the entropy profile normalization to deviate from the self-similar 
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scaling while /icm remains independent of cluster mass. If instead the similarity-breaking is 
introduced because cooling and AGN/supernovae feedback produce a mass-dependent ICM 
mass fraction, then K m&x oc T^ 3 implies /icm T^ 2 . Changing the entropy normaliza- 
tion with fixed /icm changes Lsz by ~ 10% in the M = 10 13 h~ l M Q cluster (where K max 
is increased by a factor of 2.78). Introducing a mass- dependent baryon fraction to account 
for this entropy scaling is equivalent to simply scaling the density profile and thus L sz by 
ficu/ fb- Note that this model implies an unreasonable reduction in the ICM mass fraction 
by 79% for the W z h~ l M Q cluster. 

For completeness we also considered the "threshold" cooling model from Voit et al. (2002). 
This model estimates the threshold entropy below which gas would condense by the cur- 
rent age of the universe, removes all gas below the threshold, and leaves the entropy profile 
for the remaining gas unmodified. The authors suggest that profile describes either strong 
feedback that ejects the gas below the entropy threshold to well beyond the virial radius, 
or weak feedback where the gas simply cools. Cooling of the gas above the entropy thresh- 
old is also neglected. For the M = 10 15 ^ _1 M o (K max = 2.5K l00 ) and M = lO 14 /!- 1 ^/© 
(Kmax — 1-5-Kioo) clusters, this model yielded low L sz values, but once accounted for the 
missing gas fraction, the L sz values were within 2% of the self-similar value. This reinforces 
our result that variation in the entropy profiles affect the average gas temperature only very 
weakly. The prescribed method in Voit et al. (2002) for computing the unmodified entropy 
distribution yielded a non-monotonic entropy profile for the M = 10 l3 h~ 1 M Q cluster, and 
so we did not complete the calculation for this cluster. 

In Figures 4 - 6 we compare our phenomenological density and temperature profiles with 
12 nearby, relaxed clusters from Vikhlinin et al. (2005). We use the 10 15 h~ l M & cluster val- 
ues, which are in the best agreement with the NFW fit values in Vikhlinin et al. (2005), 
cioo — 3.6 — 9.2. Our adopted value /b = 0.132 is w 25% below the universal value con- 
strained by CMB observations to 0.175 ± 0.023 (Readhead et al. 2004; Spergel et al. 2003) 
and is in reasonable agreement with the observed stellar component of clusters (Roussel 
et al. 2000; Balogh et al. 2001; Lin et al. 2003). Near t^qq, observed density and temperature 
profiles agree with the model curves with w accr = 1 (see Figures 3, 4, and 6), consistent with 
our assumptions about the outer gas, accretion pressure, and /b. In Figure 4, we also see 
the T > 5 keV observed cluster density profiles suggest an entropy profile shallower than r 11 
and/or lower ca, particularly in the core. The average temperature profile agrees well with 
our r 7 model (long dashed curve in Figure 6) outside 0.2ri o- This comparison indicates 
that our chosen parameter values cover a range of gas profiles at least as large as the scatter 
in the observations. 

Under the self-similarity assumption, p(r/r vir ) is independent of cluster mass. However, den- 
sities in the T < 5 keV clusters (see Figure 5) are significantly lower than in the hottest 
clusters (Figure 4). We devise two models to trace the lowest density profile in Figure 5. 
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The resulting temperature profiles (thick long-short dashed curve in Figure 6) agrees with 
the observed hot cluster average profile out to 0.35r„j r , and has a sharper inner peak, as 
observed in Vikhlinin et al. (2005) cool clusters. Fits 1 and 2 both have si = s 2 = 0.7 and 
ca = 7.7, the value of Cioo expected by Dolag et al. (2004) for a 10 14 /i~ 1 M Q cluster. In fit 

1 we keep the same ICM mass fraction, /icm = 0.13, as in the phenomenological models 
shown in Figures 2-4, and raise the gas entropy normalization to K max = 2.5K W0 . In 
contrast, fit 2 assumes the accretion pressure is set by gas with f h = 0.13, while the ICM 
mass fraction inside the cluster is /icm = 0.5/b, either due to ejection or increased con- 
densation. The entropy normalization in fit 2 is the self-similar value for this ICM mass 
fraction: K max = K 100 (0.5f h )~ 2 / 3 . By design, fits 1 and 2 are nearly indistinguishable in the 
X-ray observable region (r < rsoo both in density and temperature profiles). However, fit 1 
extends to 1.56r vir and contains all of the gas initially associated with the region, while fit 

2 extends only to 0.84r vir and contains half as much hot gas. This example demonstrates 
the inability of X-ray observations to constrain the distribution of non-gravitational heating 
required for similarity breaking; in fit 1, the non-gravitational heating increases the gas en- 
tropy levels and its potential energy, while in fit 2, the non-gravitational heating ejects half 
of the gas from the cluster region and leaves the remaining gas undisturbed. However, the 
SZE observables for these two models are different, and we shall return to this point in Sec. 6. 
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Fig. 2. — Phenomenological model density profiles for our 10 15 /i~ 1 M Q cluster. See Table 1 
for corresponding observables. The self-similar model is shown in each panel as the dashed 
curve. In the upper left if majf is varied between K 2 oo and SK 100 . The central density decreases 
and r ac increases as K max increases. In the upper right the entropy profile is varied. The 
central density increases with radial entropy power s 2 (outer curves), but flattens when an 
entropy core is included (inner curves). In the lower left the central density increases with 
ca for ca = 3 to ca = 11- In the lower right the accretion pressure is varied through w accr . 
The density at r ac increases with u> accr , while the central density decreases. 




Fig. 3. — Phenomenological model temperature profiles for our 1O 15 /i~ 1 M cluster scaled by 
the dark matter virial temperature Panels as in Figure 2. At fixed radius temperature 
increases with K ma , x (upper left), ca (lower left), and w accr (lower right). Temperature 
decreases with radial entropy power S2 (outer curves in the upper right), and an entropy 
core produces sharply peaked central temperature profiles (upper right). 
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Fig. 4. — All phenomenological model density profiles for our 10 15 /i _1 M Q clusters (dotted 
curves, also in Figure 2) compared to observed T > 5 keV clusters from Vikhlinin et al. 
(2005) (solid curves). The observed density profiles are consistent with w accr = 1; our 
models varying the accretion pressure (thick long-short dashed curves) do not agree the the 
observed density gradient at large radii. The observed profiles are consistent with models 
with lower ca or more entropy at small radii compared with our self-similar model. 




Fig. 5. — Observed cool clusters from Vikhlinin et al. (2005) grouped by temperature: 
2.5 keV < T < 5 keV (short dashed curves), and T < 2.5 keV (long dashed curves). The 
solid curves are two additional phenomenological models devised to trace the lowest density 
profiles in the observed region and discussed in Sec. 3.3. While the accreting gas is assumed 
to have the same baryon fraction in both models, one has 50% less gas in the ICM. These 
model density profiles are indistinguishable out to the maximum observed radius, though 
r ac is much lower in the /icm = 0.5/b model. These models demonstrate the degeneracy 
between a high entropy level and low baryon fraction in the region of gas currently X-ray 
observable. 
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Fig. 6. — Temperature profiles for our lO 15 /i _1 M model clusters (dotted curves) as in Fig. 3, 
but each normalized by their gas mass- weighted temperature within r 50Q , T 500 . The thick 
solid line is an average temperature profile for observed clusters with T > 2.5 keV from 
Vikhlinin et al. (2005). As can be seen from Fig. 3, the gradient of T/T 500 increases with 
increasing ca and decreasing K max , w accr , and S2- Outside 0.2r vir our s 2 = 0.7 model agrees 
well with the average observed profile (long dashed curve). The thick short-long dashed curve 
is the temperature profile of our two additional phenomenological models (indistinguishable 
on this plot) devised to fit low temperature cluster density profiles (see Sec. 3.3). They 
provide a good fit to the average temperature profile out to 0.35r vir . 



Finally we consider a simple polytropic model P oc p 1 in an NFW potential and fix 
7 = 1.2, in rough agreement with both simulations and observations, at least outside the 
core (see Ostriker et al. (2005), Appendix A of Voit et al. (2003), and references therein). 
The hydrostatic equilibrium equation then reduces to 

dT 7 — 1 dcpNpw 
— = -pm p 

ar 7 ar 

The NFW concentration parameters are fixed again by Eqn. 6. The free parameters in the 
model are the constant of integration in the temperature profile obtained from Eqn. 16, the 
normalization of the density profile, and the accretion radius. We assume a strong shock 
at the accretion radius and set the pressure at the accretion radius as in the models above. 
Requiring f^M{t ) gas to be contained within the accretion radius constrains the density 
normalization. However, we do not wish to impose a strong condition on the density at the 
accretion radius, since the global density of the gas accreting at t is not the density deter- 
mining the post-shock entropy in lumpy accretion; compression of gas accreting in distinct, 
dense lumps will not significantly raise the global density at r ac . We vary x ac , and the other 
free parameters are then determined by the conditions above. We exclude models where 
the post-shock density does not fall between pi and 4pi, the strong shock smooth accretion 
limit. Note that these models produce entropy profiles with L^(0) > and monotonically 
decreasing temperature profiles (see Eqn. 16). Again as shown in Appendix A we find the 
L sz signal to be extremely robust to our variation of x ac (and thereby the post-shock entropy 
normalization). 
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4. Scaling Relations 

4.1. Scalings between Observables 

While our set of models show a large variation in X-ray properties and central Compton 
parameter (see Appendix A), the total SZE luminosity is remarkably robust to the varia- 
tions introduced in our models. In Figure 7 we see that the scatter for our model clusters 
is small in ^sz/(/icM//b), indicating a tight Lsz — Ma relationship for known /icm, despite 
the consideration of a wide range of possible entropy profiles and cluster unknowns. Note 
that Lsz oc /icm makes the relation sensitive to the fraction of baryons cooled or ejected 
from the ICM. This plot corroborates the finding in McCarthy et al. (2003a) that y Q is much 
more sensitive to -preheat (or other non-gravitational physics) than is Lsz- 
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Fig. 7. — Results from our analytic baseline, preheating, and cooling models with -R'preheat 
ranging from 10 30 to 10 34 erg cm 2 g 5 / 3 , as well as all our phenomeno logical and polytropic 
models (see Sec. 3) varying the entropy profile K(r), dark matter concentration ca, accretion 
pressure, and ICM mass fraction. Crosses - 10 15 /i _1 M o cluster models, triangles - lO 14 /i _1 M 
cluster models, squares - lO 13 /i _1 M cluster models. The SZE luminosity, L S z, normalized 
by /icm/Am shows little variation within our class of models. Thus, given the ICM mass 
fraction /icm, -^sz is an excellent cluster mass indicator. The central Compton parameter, 
y , varies widely in our class of models and thus contains information on the thermal history 
of the gas. 
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Fig. 8. — Lx,cuti the X-ray luminosity outside projected radius 0.05rA, vs. central Compton 
parameter, y . Crosses - lQ 15 h~ l M & cluster models, triangles - 1O 14 /i _1 M cluster models, 
squares - 10 13 h~ 1 M & cluster models. The power law fit shown has exponent 1.2. Note that 
clusters spanning two orders of magnitude in mass and a wide variety of cluster parameters 
are well described by a single y Q — L x , C ut relation. In the inset we compare the observations 
assembled in McCarthy et al. (2003b) (crosses with error bars) to our lO 15 /i _1 M cluster 
models (triangles). The y — Lx normalization depends on mass for our models and is steeper 
than y —Lx,cut- The observed clusters are very massive, and we find good agreement between 
the observations and our lO 15 /i _1 M model clusters. 
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For the class of models studied here we also find a tight correlation between the X-ray 
luminosity outside the core, L x , C ut, and the central Compton parameter (see Figure 8). 
This is not surprising if nicM( r ) ~ 1/V 2 outside the core, so that if we ignore any radial tem- 
perature dependence, y ~ J nicyid^ and L x ~ J n 2 CM r 2 <ir are proportional to one another. 
Our scaling Lx oc y x 2 is in excellent agreement with similarly robust relations reported by 
Cavaliere & Menci (2001) and McCarthy et al. (2003a). We find the relation is closer to 
mass-independent when excluding the core X-ray luminosity, but we demonstrate agreement 
between our model and observations assembled by McCarthy et al. (2003b) in the y D — L x 
plane. 



Thus far we have focused on the effects of the gas entropy profile on ICM observables. 
In an attempt to understand the stability of our L sz signal to such wide variations, we 
compute changes in potential and thermal energies induced by changes in model parameters. 
To compute the gas potential energy, we estimate the cluster potential by an NFW profile, 
thus ignoring the deviation of the gas density profile from the dark matter profile: 



The gas thermal energy is directly proportional to L sz . Appendix A shows that changes in 
total gas energy result mostly in changes in potential energy. Moreover, fractional deviations 
of the mean energies of our models are on average larger in potential than thermal. Thus, 
even with uncertainty about the amount of energy injected into the ICM throughout its 
history, for the class of models studied here most of the injected energy manifests itself in 
potential energy. Note that our models produce total energy values differing by at least a 
factor of two for a given cluster mass. 

Consider the two model variations inducing the greatest L sz signal variation. Increasing the 
dark matter concentration parameter effectively increases the temperature associated with 
the dark matter potential, and so in that case we expect a significant change in thermal 
energy if T gas ~ Td ar k- Similarly, an increase in external pressure requires an increase in the 
internal pressure, thus driving up the thermal energy of the outer gas. 
An SZE-only observation can infer both the thermal (L S z oc M ICM 7i CM ) and potential en- 
ergy of the hot ICM. We have found a tight relation M A oc , and the total cluster mass 
is directly related to the cluster virial temperature and radius. Even allowing for the large 
variation of /icm in our models, Figure 9 shows that the projected half-luminosity radius 
of the SZE profile, scaled by L l j£ ~ r A , is well correlated with the potential energy per 



4.2. Cluster Energetics 
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particle scaled by the cluster virial temperature. Thus by measuring both L$z and rsz, one 
can estimate the level of energy injection now manifest in potential energy of the gas. 
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Fig. 9. — The potential energy per particle scaled by the cluster virial temperature Ta can 
be estimated by the SZE observable rs%L s ^ 5 . This relation holds over a mass range of 2 
orders of magnitude and over our entire range of model variations. Ta is fixed by the cluster 
mass and so is also tightly related to L sz . Thus the total potential energy per particle can 
be estimated from an SZE observation measuring r sz and L sz . Crosses - 1O 15 /i _1 M cluster 
models, triangles - 1O 14 /i _1 M cluster models, squares - 10 13 /i _1 M o cluster models. Ta,i5 
is the virial temperature and L S z, ss i5 is the self-similar SZE luminosity for a 1O 15 /i _1 M 
cluster. 



5. Implications of Models 



In this paper, we have studied the SZE properties of clusters. Motivated by observations, 
models reproducing the observed X-ray scaling relations, and the assumptions of spherical 
symmetry and hydrostatic equilibrium, we have spanned a range of plausible ICM density 
and temperature profiles. We expect that the range of parameters we have explored in our 
models encompasses the properties of real clusters. While many researchers attempting to 
understand deviations from self-similar behavior in X-ray properties of clusters have consid- 
ered a wide range of possible modifications to the gas entropy profiles, we have demonstrated 
that the SZE luminosity depends only very weakly on the shape and normalization of the 
entropy profile. While cooling can significantly reduce the L$z signal as a large gas frac- 
tion is cooled, the average temperature of the remaining gas is surprisingly robust. Thus, 
Lsz is an excellent measure of the parameter combination JicmM^ 3 . The robustness of 
L sz oc M ICM T A indicates that for the scope of possible entropy profiles and cluster parame- 
ters considered here, any significant injected energy must primarily serve to increase the gas's 
potential energy by expanding the gaseous region. However, if non-gravitational processes 
result in unaccounted cluster-mass dependent gas ejection rates or star formation efficiency, 
our L sz — M relation will be distorted. As shown by Vikhlinin et al. (2005) and discussed in 
Sec. 3, the ICM mass fraction appears lower out to r 50 o in cooler clusters. Therefore the total 
ICM mass fraction either depends on cluster mass, or the SZE profiles will be more diffuse 
in lower mass clusters. Recent observations of cluster mass/near-infrared K-band luminosity 
relation (Lin et al. 2004), estimates of intracluster light (Lin & Mohr 2004), and hydrody- 
namic simulations finding an ICM mass fraction increasing with mass (Kravtsov et al. 2005) 
all suggest the former at some level. We discuss techniques for constraining / ICM in Sec. 6. 
Of all the variations we have considered, L sz is more sensitive to the dark matter concen- 
tration parameter, the accretion pressure, and the ratio of gas to total mass in the cluster 
than the shape and normalization of the entropy profile. As shown in Sec. 3, observations 
support our understanding of dark matter concentration parameters and the gas boundary 
accretion pressure. 

In the context of these models, we have shown the most important and uncertain parameter 
of the thermal history of the gas is the ICM mass fraction, /icm; properties of the entropy 
profile induce variations in Lsz clearly bounded by 10% while Lsz oc /icm- Lsz shows 
variations of « 8% from the expected intrinsic scatter in ca (Dolag et al. 2004). Knowledge 
of ca and its scatter can be acquired by thorough N-body simulations, and the analysis 
of Vikhlinin et al. (2005) finds agreement between ACDM concentration parameters and 
observations of nearby, relaxed X-ray clusters. Finally, large deviations of the bounding 
pressure from our expectations could potentially significantly alter Lsz- The comparison to 
observations discussed in Sec. 3 indicates at least mild agreement with our expected values. 
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Thus we expect the scatter induced by effects discussed here to be bounded by 10%. We 
caution, however, that our assumption of hydrostatic equilibrium may not hold in the out- 
skirts of the cluster (Voit et al. 2002; Thomas et al. 2002). Furthermore, our modeling is far 
from extensive: we have ignored cluster asphericity, thermal conduction, intracluster mag- 
netic fields, turbulent support, the presence of a relativistic fluid, and important dynamical 
events such as large mergers and quasar blast waves. In an analytic gas model, Ostriker 
et al. (2005) demonstrate that triaxiality and substructure in simulated dark matter halos 
do induce some scatter in the cluster observables; they find cr hiyi0Q m 0.3 for the y wo — M 1Q0 
relation when M 100 > 10 14 M Q . Motl et al. (2005) find in hydrodynamic simulations includ- 
ing star formation and supernovae feedback that while major mergers may increase y Q by 
up to a factor of 20, they do not drastically increase y^oo (see also Randall et al. (2002)). 
Furthermore, the level of scatter is similar to our findings: 80% of the simulated cluster mass 
estimates from the 1/500 — M500 relation lie within +15% to —10% of the true cluster mass. 
The scatter level is nearly constant with redshift back to at least z = 1.5, supporting our 
optimistic view that merging has a small effect on Lsz- Thus we are optimistic that Lsz will 
be a useful indicator of cluster mass, while SZE profiles coupled with measurements in other 
wavebands can help us constrain the magnitudes and sources of non-gravitational physics 
important to understanding clusters. 

Our models have indicated a number of tight relations between cluster properties and ob- 
servables: 

. ^^i.exio-Mpc 2 ^,^^ 573 

• Lx,cut oc y l Q 2 with a normalization independent of cluster mass 

• The net energy injected into the bound ICM can be estimated from the observable 

— 1/5 

r szL sz as in Figure 9 and discussed in Sec. 4.2. 

These relations should now be examined in hydrodynamic simulations, where additional 
physics can be modeled and departures from hydrostatic equilibrium and spherical symme- 
try can be examined. 



6. Constraining /i CM 

Within both our models and current observations, we have shown the most important 
uncertain parameter in the Lsz — Ma relation is the ICM mass fraction, /icm- Can we 
observationally determine /icm as a function of both cluster mass and redshift? 
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SZE surveys can potentially break the /icm — Ma degeneracy inherent in Lsz measurements 
through the detection of the kinetic Sunyaev-Zel'dovich effect (kSZ). A cluster's peculiar 
velocity v produces a CMB temperature change given by (Sunyaev & Zeldovich 1972) 

' ' Z = (^Thomson f dln e (l) ( ^ - ) . (18) 



T, 



o 



The kSZ temperature change averaged over cluster i, 5i/T Q = —o 'Thomson^ T e,i ^Los.i/c, is pro- 
portional to the total number of cluster electrons and the cluster peculiar velocity projected 
along the line of sight. Hernandez-Monteagudo et al. (2005) discuss an unbiased estimator 
of Sf 

Si = 5TcMB,i + STts Zti + N + STlf z . + T kSZ ,. (19) 

The variance of Si will therefore have contributions from the desired kSZ signal and a noise 
term from the residual CMB and thermal SZ signals after subtraction, 8T^ e £ IB i + 5T t sz,i, the 
instrumental noise, N, and internal cluster motions, ST^ Zi . 



<(!)>= ^sonNUi^f + °L- (20) 

With a good understanding of of es and given the concordance AC DM prediction for (v^ os ), 
one may deduce the total electron content and thereby Micm for a sample of clusters in a 
certain Lgz and redshift range, thus constraining ficwiM, z). This technique will be par- 
ticularly applicable to SZE surveys, where ST^z and Lsz can be measured within the same 
projected area to constrain /icm without extrapolation in radius. We caution, however, that 
the kSZ effect has not yet been detected, and so the systematic errors associated with this 
measurement are uncertain. 

As we have demonstrated in Sec. 3.3, two different models for the cooling and feedback en- 
ergy distribution can reproduce the density and temperature profiles of observed cool clusters 
deviating from self-similarity. Though they differ in /icm by a factor of 2, these models are 
indistinguishable out to ps 0.85r V i r , well beyond the typical X-ray detectable region inside 
~ 0.5r v i r . In the first (fit 1), non-gravitational heating is well distributed and raises the 
entropy level of all the cluster gas, and thus increases its potential energy (see discussion in 
Sec. 4.2). In the second (fit 2), half of the ICM is either cooled or ejected from the clus- 
ter, while the remaining hot gas maintains the properties of a gravitationally heated ICM. 
These two cluster models do have distinct SZE observables: fit 1 has L S z = 0-88L S z, ss and 
r sz = 0.49, while fit 2 has L sz = 0.57L 5ZjSS = 1.13L 5z , ss /(/ ICM // b ) and r sz = 0.32. Both 
models cause < 13% deviation of L S z / {hcu/ fb) from the self-similar value. We emphasize, 
however, that relating the observed thermal SZE profile to /icm would still require an as- 
sumption about the behavior of p{r) or T(r) in the outer regions of the cluster. Since X-ray 
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observations are not sensitive to the outer regions of the cluster, alone they can only place 
a lower limit on the number of baryons in the ICM and cannot distinguish between these 
two heating models. Ma can be estimated using other means, such as gravitational lensing 
or assuming the X-ray observable gas is in hydrostatic equilibrium. These follow-up mea- 
surements could also be used to calibrate the L$z — Ma relation and as a consistency check 
with the kSZ results. Because systematics will vary with the technique, multi-wavelength 
cluster observations will enhance our ability to minimize our assumptions and improve our 
accuracy in the calibration of the L S z — Ma relation, as well as bring new understanding to 
the physics of cluster gas. 



7. Conclusions 

SZE surveys will soon identify thousands of clusters and measure their photometric 
redshifts, so we will be able to estimate the gas potential and thermal energy with L sz 
and r S z, as well as the cluster mass (see Fig. 7 and Fig. 9). Detection of the kSZ signal 
in these SZE surveys will also constrain the ICM mass fraction, /icm(M, z). Follow-up ob- 
servations in X-ray, optical, and radio can provide checks on the survey Lgz — M relation 
using independent mass determinations, and could further constrain ICM models through 
AGN activity-gas energy correlations. Coupling sufficiently high resolution SZ and X-ray 
observations could over-constrain the gas density and temperature profiles, and thus provide 
a consistency check of deprojection techniques. Sensitive SZE profiles should extend farther 
into the outer regions of the cluster than X-ray observations allow. 

This paper has focused on SZE observables from z = clusters, though SZE surveys will 
find most of their clusters at z ~ 0.6. However, the models studied here support the self- 
similar prediction T gas ~ T^ark even under considerable variations of the total gas energy. 

1 /3 

We expect T^ark to scale with the overdensity (A{z)p cr (z)) 1 , and the dark matter concen- 
tration parameter will decrease with redshift as suggested by Eqn. 6: by a factor of ~ 1.5 at 
z — 1. Appendix A suggests that changes in L^z induced by a decrease in ca with redshift 
is limited to less than 15%. The range of entropy normalizations examined in this paper is 
larger than the expected characteristic entropy decrease by ~ 20% at z = 1 under self-similar 
evolution. Because we have found a very weak dependence of Lsz on the entropy profile, 
we are optimistic that the redshift evolution of the Lg Z — M relation will be determined 
by the more easily understood dark matter evolution; this assumption has been verified in 
numerical simulations including galaxy formation (Nagai 2005). In this case the Lsz — M 
evolution is characterized by relatively few parameters, and so can hopefully be internally 
calibrated in SZE surveys (Hu 2003; Majumdar & Mohr 2004; Lima & Hu 2004, 2005; but 
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see also Francis et al. 2005). However, we caution that baryon physics may induce a redshift 
dependent ICM mass fraction ficM.(M,z) that produces a deviation from simple self-similar 
redshift evolution. The observed ICM mass fraction in clusters appears to depend on cluster 
mass and radius. In Sec. 6 we proposed the use of the kSZ signal to constrain ficM.(M,z). 
Understanding the variation in ICM mass fraction will deepen our understanding of cluster 
physics and enable the use of clusters as cosmological probes. Finally, we may hope to probe 
the gas's thermal history by measuring the thermal and potential energies as a function of 
redshift through SZE observation. 
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A. Appendix A 

We list SZE and X-ray observables computed from the models considered in the paper. 
r ac is the accretion shock radius, given in units of the dark matter virial radius, r&- ">"sz 
and r x are the radii projected on the sky enclosing half the SZE/X-ray luminosity, and also 
given in units of r&. y is the unitless Compton parameter measured along the line of sight 
through the center of the cluster. We report Lsz in terms of Lsz, ss , the value expected from 
the self-similar collapse model: 

T ^Thomson fbM & , n -4,, 2 / A 1 \ 

where /b is the universal baryon fraction and Ta is the characteristic temperature of the 
dark matter potential. Lx is the total X-ray luminosity (see Eqn. 13), and Lx, C ut is com- 
puted by integrating the X-ray luminosity outside a projected radius of 0.05rA- The X-ray 
luminosities were computed with the cooling function found in Sutherland & Dopita (1993), 
which includes emission at all wavelengths. We also compute the average potential (PE), 
thermal (3/2/cT), and total energy (E) per particle in keV according to the approximation 
made in Section 4.2. 1 — /icm/ /b is the gas fraction that has been cooled or ejected from 
the cluster. 

The tables below list first our baseline model without heating and cooling, followed by 
models including preheating and cooling with various values of preheat given in units of 
erg cm 2 g~ 5 / 3 . The parameters for the "self-similar" phenomenological model for the 1O 15 /i _1 M 
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cluster are ca = 8.5, u> accr = 1 (the "fudge factor" on the mass accretion rate), s± = s 2 = 11 
(entropy exponents), x core = (entropy core radius), K max = K 100 , K min = 0.01K max , and 
/b = 0.13. The "self-similar" parameters are the same for the 1O 14 /i _1 M and 10 13 h~ 1 M Q 
clusters except ca = 10.4 and ca = 12.6 as used in the preheat model of Voit et al. (2003). 
Fit 1 in Table 1 (see Sec. 3.3 for discussion) has parameters ca = 7.7, u> accr = 1, s± — S2 — 0.7, 
Xcore = 0, fT max = 2.5Kioo, K min = 0.01K max , / b = 0.13. Fit 2 has parameters c A = 7.7, 

W accr = /b//lCM, Si = S 2 = 0.7, X core = 0, K max = i^loo(/lCM//b) _2/3 , Kmin = 0.01K max , 

/ b = 0.13, and h CM / h = 0.5. 
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Table I. l^h^M,, Cluster Models 



Model 




»*sz 


Vo 


Lsz 


r x 


Lx 


Lx,cut 


PE 


ffcT 


E 


1 - /icm/A. 




(•■a) 


(Ta) 


(xlO 4 ) 


(Lsz,ss) 


(ta) 


(10 44 


erg s _1 ) 




(keV) 






no heat /cool 


0.893 


0.255 


3.6 


0.963 


0.0293 


55. 


18.1 


-29.5 


9.33 


-20.2 





^"preheat = 10 34 


1. 


0.297 


1.28 


0.97 


0.123 


11.5 


9.84 


-25.5 


9.4 


-16.1 





K ph = 5 x 10 33 


0.946 


0.272 


1.94 


0.965 


0.0769 


21.1 


14.8 


-27.5 


9.36 


-18.1 





K ph = 4 x 10 33 


0.935 


0.268 


2.23 


0.963 


0.0635 


26. 


16. 


-28. 


9.34 


-18.7 





K ph = 3.5 x 10 33 


0.93 


0.265 


2.45 


0.961 


0.0549 


29.9 


16.3 


-28.3 


9.32 


-19. 





K ph = 3 x 10 33 


0.924 


0.263 


2.81 


0.957 


0.0428 


36.8 


16.1 


-28.6 


9.29 


-19.3 





K ph = 2 x 10 33 


0.913 


0.263 


2.65 


0.939 


0.0476 


31.4 


15.1 


-28.5 


9.34 


-19.1 


0.0252 


K ph = 10 33 


0.902 


0.261 


2.66 


0.928 


0.0468 


31.6 


15.1 


-28.6 


9.36 


-19.3 


0.0378 


K ph = 5 x 10 32 


0.896 


0.259 


2.69 


0.923 


0.046 


32. 


15. 


-28.7 


9.37 


-19.4 


0.0442 


K ph = 10 32 


0.892 


0.258 


2.71 


0.919 


0.0452 


32.4 


15. 


-28.8 


9.37 


-19.5 


0.049 


K ph = 10 31 


0.891 


0.258 


2.67 


0.918 


0.0464 


31.8 


15. 


-28.8 


9.37 


-19.4 


0.0501 


K ph = 10 30 


0.891 


0.258 


2.67 


0.918 


0.0465 


31.7 


15. 


-28.8 


9.37 


-19.4 


0.0502 


self similar 


0.806 


0.227 


3.41 


0.97 


0.0437 


51.4 


23.2 


-31.2 


9.46 


-21.7 





c A = 3 


0.887 


0.312 


0.87 


0.762 


0.092 


14.5 


10.2 


-24.1 


7.32 


-16.8 





c A = 7 


0.819 


0.243 


2.54 


0.923 


0.0509 


38.6 


19.5 


-29.3 


8.89 


-20.4 





C A = 11 


0.781 


0.205 


5.08 


1.04 


0.0356 


74.6 


28.3 


-33.4 


10.1 


-23.2 





tUaccr = 0-28 


0.906 


0.192 


3.82 


0.855 


0.0406 


64.9 


27.5 


-32. 


8.32 


-23.6 





^accr — 3.5 


0.65 


0.268 


3.03 


1.39 


0.0583 


40.2 


22. 


-30.8 


13.3 


-17.5 





si = S2 = 0.7 


0.944 


0.295 


1.5 


0.968 


0.105 


12.5 


9.41 


-25.8 


9.42 


-16.4 





si = s 2 = 1.5 


0.669 


0.167 


6.42 


0.946 


0.0307 


195. 


53.3 


-37. 


9.14 


-27.9 





si = 0, x c = 0.1 


0.813 


0.231 


2.17 


0.982 


0.0857 


28.3 


21.3 


-29.7 


9.45 


-20.3 





X min = 3.5 x 10 33 


0.813 


0.23 


2.34 


0.979 


0.0783 


31.4 


22.5 


-30.2 


9.5 


-20.7 





^"max = ^200 


0.719 


0.204 


3.93 


0.97 


0.0419 


67.7 


29.6 


-33. 


9.48 


-23.5 





A"max = 1.5A"lO() 


0.988 


0.274 


2.64 


0.964 


0.047 


31.2 


14.9 


-27.8 


9.31 


-18.5 





Kmax = 2iiTl00 


1.15 


0.316 


2.21 


0.951 


0.0493 


22.3 


11. 


-25.7 


9.24 


-16.5 





ffmax = 2.5KlO() 


1.3 


0.353 


1.93 


0.936 


0.0511 


17.2 


8.73 


-24.2 


9.16 


-15. 





^max = 3-fSTlOO 


1.43 


0.385 


1.71 


0.926 


0.0526 


13.5 


7. 


-22.5 


8.92 


-13.6 





Cooling threshold 


0.993 


0.265 


1.78 


0.871 


0.0829 


18.2 


13.2 


-28.3 


9.68 


-18.6 


0.128 


Sec. 3.3 fit 1 


1.56 


0.488 


0.65 


0.882 


0.144 


3.01 


2.45 


-18.3 


8.62 


-9.67 





Sec. 3.3 fit 2 


0.838 


0.32 


0.66 


0.565 


0.130 


2.91 


2.33 


-25.6 


10.9 


-14.7 


0.5 


7 = 1.2, p 2 = 3.96pi 


0.875 


0.25 


1.92 


0.981 


0.0927 


21.5 


16.1 


-28.5 


9.51 


-19. 





7 = 1.2, P 2 = 2.53pi 


1.2 


0.341 


1.23 


0.953 


0.108 


9.25 


7.31 


-23.7 


9.24 


-14.4 





7 = 1.2, p 2 = 1.74pi 


1.5 


0.421 


0.914 


0.917 


0.117 


5.26 


4.26 


-20.5 


8.89 


-11.6 





7 = 1.2, p 2 = 1.05pi 


1.83 


0.483 


0.758 


0.866 


0.119 


3.64 


2.96 


-18.3 


8.4 


-9.86 
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Table 2. 10 m h~ l M & Cluster Models 



Model 


r ac r S z 




r x 


Lx L x ,cut PE ffcT E 


1 - /lCM//b 




(i-a) Oa) 


(xlO 5 ) (L SZ , SS ) 




(10 43 ergs- 1 ) (kcV) 





no heat /cool 



A- 



preheat 



10 



34 



ph 



A ph 

A-ph 
A ph 
A ph 
A ph 
A ph 
Ap h 
Aph 
ph 



5 x 10 33 
: 4 x 10 33 
: 3.5 x 10 33 
: 3 x 10 33 

: 2 X 10 33 

: 10 33 

: 5 X 10 32 

: 10 3 

: 10 : 



32 
31 



A, 



10 30 
self similar 
c A = 4.5 
c A = 8.75 
ca = 13 

Waccr = 0.28 
W^accr — 3.5 

si = S2 = 0.7 
si = S2 = 1.5 
si = 0, x c = 0.1 
A min = 3.5 x 10 33 
A max = A200 
A max = 1.5Aioo 
A max = 2Kioo 
A max = 2.5Aioo 
A max = 3Aioo 

Waccr = 1.78 
J1 max 
A ma x ^ 

A- 

Cooling threshold 
7 = 1.2, p 2 = 3.87pi 
7 = 1.2, p 2 = 2.64pi 
7 = 1.2, p 2 = 1.74/31 
7 = 1.2, p 2 = 1.03pi 



, 7^ 
1 A 

T 2/3 
A ' 
T l/2 



1.03 

1.47 

1.26 

1.22 

1.19 

1.17 

1.11 

1.05 

1.03 

1. 

1. 

1. 

0.838 

0.919 

0.85 

0.813 

0.912 

0.713 

1. 

0.681 

0.85 

0.95 

0.744 

1.04 

1.21 

1.38 

1.54 

0.787 

1.1 

0.838 

0.97 

0.975 

1.25 

1.55 

1.83 



0.227 

0.436 

0.325 

0.3 

0.288 

0.276 

0.251 

0.249 

0.244 

0.241 

0.24 

0.24 

0.192 

0.261 

0.205 

0.173 

0.168 

0.236 

0.263 

0.133 

0.198 

0.261 

0.173 

0.229 

0.26 

0.287 

0.31 

0.21 

0.24 

0.192 

0.251 
0.214 
0.263 
0.309 
0.332 



5.18 

0.446 

0.888 

1.1 

1.24 

1.45 

2.28 

2.95 

3.05 

2.95 

2.93 

2.99 

4.91 

1.47 

3.74 

7.07 

5.39 

4.26 

1.98 

9.4 

2.91 

1.21 

5.65 

3.83 

3.21 

2.81 

2.54 

4.62 

3.6 

2.28 

1.34 

2.87 

2.13 

1.7 

1.53 



0.932 

0.846 

0.91 

0.92 

0.924 

0.929 

0.932 

0.852 

0.828 

0.808 

0.804 

0.804 

0.957 

0.753 

0.909 

1.03 

0.896 

1.19 

0.936 

0.947 

0.968 

0.948 

0.964 

0.936 

0.916 

0.896 

0.874 

1.03 

0.929 

0.444 

0.705 
0.958 
0.923 
0.883 
0.845 



0.0227 

0.326 

0.197 

0.165 

0.148 

0.129 

0.0821 

0.0439 

0.0408 

0.0425 

0.0429 

0.041 

0.0357 

0.0617 

0.04 

0.0303 

0.0346 

0.0401 

0.0929 

0.0287 

0.0809 

0.186 

0.034 

0.0387 

0.0409 

0.0428 

0.0444 

0.0369 

0.0395 

0.0357 

0.111 

0.0751 

0.0821 

0.0866 

0.0871 



54.6 

I. 29 
3.48 
4.72 
5.66 
7.04 
13.8 
17.5 
18.4 
17.4 
17.2 
17.8 
48.4 
18.4 
38. 
66.3 
58.9 
35.5 
9.41 
200. 
23. 
6.35 
63.9 
29.8 
21.1 
16.3 
13.5 
42.5 
26.3 
10.4 

5.42 
20.1 

II. 3 
7.3 
5.95 



13.2 

1.24 

3.23 

4.26 

5.02 

6.07 

10.1 

7.95 

7.9 

7.7 

7.66 

7.67 

18.4 

10.5 

16. 

21.4 

21.5 

15.2 

6.67 

45.7 

16.8 

5.85 

23.3 

12.2 

9.02 

7.21 

6.14 

16.8 

10.9 

3.97 

4.43 
13.5 
7.97 
5.26 
4.3 



-6.63 

-3.9 

-4.87 

-5.15 

-5.31 

-5.5 

-5.98 

-6.13 

-6.23 

-6.28 

-6.29 

-6.3 

-7.18 

-5.75 

-6.76 

-7.65 

-7.35 

-6.97 

-5.75 

-8.75 

-6.83 

-5.67 

-7.62 

-6.43 

-5.91 

-5.54 

-5.29 

-7.08 

-6.25 

-7.18 

-6. 

-6.56 

-5.75 

-5.12 

-4.78 



1.94 

1.77 

1.9 

1.92 

1.93 

1.94 

1.95 

1.96 

1.97 

1.97 

1.97 

1.97 

2. 

1.58 

1.88 

2.14 

1.86 

2.46 

1.94 

1.98 

2.02 

1.99 

2.03 

1.95 

1.9 

1.86 

1.84 

2.15 

1.94 

2. 

2.05 
2. 

1.93 
1.84 
1.77 



-4.68 

-2.13 

-2.97 

-3.23 

-3.38 

-3.55 

-4.03 

-4.18 

-4.26 

-4.3 

-4.31 

-4.33 

-5.17 

-4.17 

-4.88 

-5.5 

-5.49 

-4.5 

-3.8 

-6.77 

-4.81 

-3.69 

-5.59 

-4.48 

-4.01 

-3.68 

-3.45 

-4.93 

-4.32 

-5.17 

-3.95 
-4.56 
-3.82 
-3.27 
-3.01 











0.0904 

0.12 

0.144 

0.149 

0.15 



































0.536 

0.282 
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Table 3. lO^h^M^ Cluster Models 



Model 


f ac 


rsz 


Vo 


Lsz 


r x 


Lx 


Lx,cut 


PE 


ffcT 


E 


1 - /lCM//b 




(r-A) 




(xlO 6 ) 


(Lsz,ss) 




(10 42 


erg s _1 ) 




(keV) 






I1U IlCeLl j l_-U(Jl 


1.15 


0.196 


7.58 


0.942 


0.0194 


94.4 


17.8 


-1.52 


0.424 


-1.1 


o 


k , — ^ v 1 n 33 

^preheat — X 1U 


2 43 


u.uuu 


n 9 


u.uuo 


n 7zl^ 






n 

-u.ooy 


O 9QR 


-0 241 


n 


K . - 4 y 1 n 33 


2.06 


0.574 


0.273 


0.717 


0.622 


0.717 


0.706 


-0.632 


0.323 


-0.309 


o 


if — o c v 1 O 33 


1.99 


0.517 


0.324 


0.717 


0.554 


0.901 


0.885 


-0.679 


0.323 


-0.357 


o 


^ ■ - 3 v in 33 


1.79 


0.464 


0.404 


0.75 


0.483 


1.22 


1.2 


-0.751 


0.337 


-0.413 


o 


K ph = 2 x 10 33 


1.43 


0.347 


0.764 


0.853 


0.312 


2.91 


2.78 


-0.967 


0.384 


-0.584 





K ph = 10 33 


1.21 


0.214 


4.85 


0.769 


0.0318 


35.5 


11.9 


-1.44 


0.422 


-1.02 


0.18 


X ph = 5 x 10 32 


1.12 


0.238 


2.69 


0.642 


0.0387 


10.6 


4.51 


-1.34 


0.429 


-0.907 


0.326 


X ph = 10 32 


1.02 


0.223 


2.53 


0.579 


0.0424 


9.51 


4.25 


-1.4 


0.437 


-0.96 


0.403 


K ph = 10 31 


1. 


0.22 


2.58 


0.566 


0.0394 


9.82 


4.17 


-1.41 


0.439 


-0.974 


0.419 


X ph = 10 30 


1. 


0.22 


2.55 


0.565 


0.0404 


9.62 


4.16 


-1.41 


0.439 


-0.974 


0.42 


self similar 


0.856 


0.163 


7.13 


0.985 


0.0303 


77.6 


25.6 


-1.67 


0.444 


-1.22 





c A = 7 


0.919 


0.21 


2.82 


0.809 


0.0438 


42.3 


19.3 


-1.4 


0.364 


-1.03 





c A = 11.5 


0.863 


0.17 


6.07 


0.954 


0.0322 


68.9 


24.2 


-1.6 


0.425 


-1.18 





ca = 16 


0.831 


0.146 


10.6 


1.07 


0.026 


102. 


28.6 


-1.8 


0.484 


-1.32 





Waccr = 0.28 


0.912 


0.148 


7.68 


0.956 


0.0301 


90.9 


29.4 


-1.71 


0.427 


-1.28 





^accr — 3.5 


0.762 


0.197 


6.31 


1.1 


0.0313 


59.1 


20.6 


-1.62 


0.501 


-1.12 





si = s 2 = 0.7 


1.05 


0.236 


2.69 


0.941 


0.0906 


12.9 


8.9 


-1.31 


0.425 


-0.885 





si = S2 = 1.5 


0.681 


0.109 


13.7 


0.993 


0.0276 


334. 


66.8 


-2.07 


0.446 


-1.62 





si = 0, x c = 0.1 


0.875 


0.172 


3.91 


0.995 


0.0809 


32.1 


23.3 


-1.58 


0.449 


-1.13 





^max = ^200 


0.756 


0.149 


8.18 


0.998 


0.0287 


102. 


32.1 


-1.77 


0.45 


-1.32 





ffmax = 1.5A"l0() 


1.08 


0.193 


5.63 


0.954 


0.0335 


48.4 


17.5 


-1.51 


0.432 


-1.08 





Kmax = 2i^l00 


1.26 


0.217 


4.72 


0.927 


0.0361 


34.2 


13.1 


-1.39 


0.417 


-0.971 





^"max = 2.5A"l0() 


1.44 


0.237 


4.13 


0.902 


0.0383 


26.3 


10.6 


-1.3 


0.406 


-0.895 





Kmax = SKlOO 


1.6 


0.253 


3.69 


0.881 


0.0402 


21.2 


8.87 


-1.23 


0.393 


-0.833 





Waccr = 3.16 


0.769 


0.193 


6.35 


1.09 


0.0311 


60. 


20.7 


-1.61 


0.489 


-1.12 





K r*j T 2 / 3 

I ^ max J A 


1.54 


0.247 


3.89 


0.889 


0.0394 


23.5 


9.68 


-1.27 


0.402 


-0.865 





ft" ~ t-" 2 / 3 

^inax A 5 


0.856 


0.163 


1.54 


0.212 


0.0303 


3.6 


1.19 


-1.67 


0.444 


-1.22 


0.785 


























7 = 1.2, p 2 = 3.87pi 


1.05 


0.177 


4.49 


0.981 


0.0643 


35. 


21.3 


-1.55 


0.442 


-1.11 





7 = 1.2, p 2 = 2.72pi 


1.3 


0.205 


3.61 


0.946 


0.0684 


22.9 


14.4 


-1.4 


0.426 


-0.973 





7 = 1.2, p 2 = 1.80pi 


1.58 


0.23 


3.06 


0.909 


0.0711 


16.6 


10.7 


-1.28 


0.409 


-0.873 





7 = 1.2, P2 = 1.06pi 


1.83 


0.24 


2.85 


0.881 


0.0715 


14.4 


9.32 


-1.22 


0.397 


-0.827 
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